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We consider the problem of temperature chaos in mean-field spin-glass models defined on random 
lattices with finite connectivity. By means of an expansion in the order parameter we show that 
these models display a much stronger chaos effect than the fully connected Sherrington-Kirkpatrick 
model with the exception of the Bethe lattice with bimodal distribution of the couplings. 
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^ ■ I. INTRODUCTION 

(N '■ 

O , Chaos in temperature is very old problem in spin-glass theory and has received a lot of attention over the years in 
1^ ■ connection with many different problems. It has been studied by various approaches including: i) scaling arguments 
Ph j and real space renormalization group analysis [1], analytical and numerical studies on mean-field models [2-6] (see 
. also [7] for a recent review), analytical and numerical studies on the finite-dimensional Edwards- Anderson model 
[8-10], analytical and numerical studies on elastic manifolds in random media [11]. The problem of the temperature 
dependence of the Gibbs measure of the Random Energy Model has been also investigated both in the physics 
literature [18] and in the mathematical physics literature [19]. It is also believed that the chaos picture is suitable 
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I ] to understand the surprising rejuvenation and memory effects observed in the dynamics of real spin-glasses, see for 
^ ' example [12-17] and references therein. Furthermore, a detailed understanding of this problem would probably shed 
^ light on the success of the Parallel Tempering procedure, which is nowadays considered an essential ingredient to 
(_i • achieve thermalization in numerical spin-glass simulations [21, 22]. 

^ ' In this paper we consider the problem of chaos in temperature in various mean-field spin-glass models with finite 
, connectivity. We show that these models display a much more pronounced chaos effect with respect to the fully 
connected Sherrington-Kirkpatrick (SK) model, with the notable exception of Bethe lattices with bimodal distribution 
of the couplings; moreover we will analyze the dependence of chaos on the various parameters. Chaos is much more 
pronounced if the system is locally heterogenous. It is possible that some of these results also hold in finite dimensional 
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Q ' models. 



II. SOME DEFINITIONS 



> 

■ In the problem of temperature chaos one is interested in the correlations between the thermodynamically relevant 
' configurations at different temperatures for a given a general spin-glass model defined by a random Hamiltonian of 
N spins Hj{a}. In particular one would like to evaluate the probability Pj^^'^{q) of observing an overlap q if we 
' extract two configurations according to their Boltzmann weights from systems with the same quenched Hamiltonian 
, but different temperatures: 

ff.fe _ EMM'^(^g-E.^.r,)exp[-/3ig,,M-/32gj{r}] 

' j:Mi.}^M^l3iHj{a}-l3,Hj{r}] ^> 

$H ' where Pi =T^ ^ and the overline represents average with respect to different Hamiltonians i7j{r}. 
. 5^ , A peculiar feature of Replica-Symmetry-Breaking theory [20] is that if the two systems have the same temperature 
the function Pj^^^{q) in the low temperature phase has a support between —qsA and qsA, where qsA is the so-called 
Edwards- Anderson parameter (for a two spin interaction in absence of magnetic field). In particular the disorder 
average P{q) — P.j{q) is given by P{q) = dx/dq where q{x) is a continuous function between zero and qsA- 

The problem of chaos in temperature concerns the function Pj^^^[q), in particular we say that there is chaos if 

P^j''Hq) = 5{q) (2) 

i.e. if Pj^^^{q) has a support concentrated on g = and that there is no chaos otherwise. In the nutshell, if chaos is 
present, the equilibrium configurations at one temperature are quite different from those at a different temperature. 
It is clear that chaos may have a dramatic effect of the dynamics after a temperature shift. 

The problem of chaos is intrinsically related to the disordered nature of these systems, being trivial in non-disordered 
models like a ferromagnet that have a translational invariant order parameter. The problem of chaos in temperature 
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in the various mean-field models has been investigated intensively over the years. Today we know that there are 
fuU-RSB models that do not display chaos in temperature [4] and fuU-RSB models that do have chaos including 
notably the SK model [5]. Similarly there are IRSB models that do have chaos in temperature and models that do 
not [6]. 

If chaos is present, we would like to quantify it, also to understand the finite volume effects (or finite time effects 
in the dynamics). More precisely we would like to know the free energy increase that happens if we constrain one 
system near to the other (or more generally at an overlap q). This free energy increase can (as usually) written also 
as large deviations function for the distribution of the overlap between systems at different temperatures: it can be 

computed by studying two coupled systems. 

The free energy of two systems forced to stay at a fixed overlap q is given by: 



i^i2(<Z,/3i,/32) = -^ln slNq-J2'^^r^eM-PlHJ{a}-/32HJ{T}]. 

{T}{a} \ i J 

The function Fi2(g, ,62) must be larger than or equal to the free energies of the two unconstrained systems and 

the relevant quantity is the free energy shift AFi2{q, /32) = -^12(9, ^^1,^2) ~ — F{02)- Indeed if this quantity 

is greater than zero it follows that the large deviations of the overlap are given by: 

P^'^^q) c< exp[-iVA^^i2(g,/3i,^2)] • (3) 

If the function P^^^^^(q) has support on some non-zero values of q then the free energy shift must vanish 
^-^12(9, /?!, /32) = 0. The opposite in general is not true, i.e. a vanishing free energy difference does not neces- 
sarily imply a non-zero P^j^^^{q) as was unexpectedly discovered in the case /3i = (32 for the spherical SK model 
[7, 30]. ' 

In the following we will consider the constrained free energy functional Fi2{q, Pi, (32) averaged over the disorder 
because it is usually assumed that this quantity, (and correspondingly the large deviations) does not fiuctuate in the 
large N limit. We could also define Pfliq, /3i, /32) the free energy of the replica a if we constrain the replica a to stay 
at a fixed overlap from replica r, the replica t being at equilibrium [27]: 



^Q, « lj:{r}^M-p2Hj{T}]lnj:{.}SiNq-j:^aiTi)eM-^iHj{a}] 

^12{1j Ply ""2) — — ^ f — a IT ( ^1 • 

N T,{t}^M-I32Hj{t}] 

In the first definition everything was symmetric in the two replicas and forcing the system to have a non zero 
overlap we push out of equilibrium both replicas. On the contrary in the second definition we look to the probability 
of a, when it is constrained to stay at a fixed overlap with r is a quenched configuration at equilibrium. It is evident 
that convexity implies that 

Fi2iq,Pi,p2)<F^^{q,pi,p2) , 

so that if Fi2{q, (3i, (32) displays chaos, chaos is present also in F^{q, (3i, (i2)- It is interesting that we can obtain the 
internal energy as function of q by performing a derivative with respect to 02- In presence of chaos the quantity 

AE{pu(32) = E%{qEA, 01,(32) - Fi1(0, A, ^2) , (4) 

should have the meaning of the energy that is slowly released after a sudden quench from a high temperature. 
Unfortunately these slow relaxations are too small to be observed experimentally. 

In this paper we compute Fi2{q, (3i, P2), the computation of F^{q, Pi, (32) could be cone in a similar way: the two 
functions display a similar qualitative behaviour. 

In order to compute the free energy shift AFi2{q, (3i, (32) is is convenient to consider the following coupled free 
energy: 



1 

'N' 







N 


In exp 


-(3iHj{a} - P2Hj{t} H 








i=l 



(5) 



where /3j is the inverse temperature of the corresponding system and A'' is the total number of spins in the system. 

When the coupling term e vanishes the above quantity is simply the sum of the averages free-energies (times /3) at 
inverse temperatures (3i and (32 , therefore the relevant physical information concerning chaos is given by the difference 
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AFi2(e, /32) = Fi2{e, di, 02) - Fi{0i) - ^2(/32)- In the thermodynamic hmit the two functions Fi2{q, pi, P2) and 
A2(ei/3i,/32) becomes Legendre transform of each other through the following relationships: 



Fi2{e) = Fi2{q*) - qe , 

and 

Fi2iq) = Fi2ie*) + qe , 



dF 
dq 



dF 
Ik 



(6) 

q—q" 



q- (7) 



III. CHAOS IN TEMPERATURE IN THE GENERALIZED SK MODEL 

In this section we consider the problem of chaos in temperature in the context of the generalized SK model, that 
is a spin-glass model whose free energy is given by the extremization of a functional of an n x n matrix Qab with 
generic coefficients. Over the years many different spin-glass models have been mapped over a generalized Sherrington- 
Kirkpatrick model e.g. the Edwards- Anderson model in an expansion in large dimension [28] and at fixed dimension 
in the loop expansion above the upper critical dimension £)„ = 6 [29]. In particular such a mapping was used by 
Kondor to assess chaos in temperature in the Edwards- Anderson model [9, 10]. In the next section we will use it to 
study spin-glass models on the Bcthc lattice. 

In all the aforementioned spin-glass models Fi2{e, Pi, P2) can be computed in the replica framework (see e.g. [7]), 
i.e. considering a set of n replicas of the two coupled systems and then sending n to zero. By means of standard 

manipulations one obtains the free energy as a functional Fi2{Q, e) over a 2n x 2n matrix Q = ^ ) ' ^^^'^^ 

QijQi and P are nx n matrices. The constrained free energy is obtained extremizing the functional with respect 
to the order parameter Q at a given value of e. The functional may be very complicated but one can get a more 
tractable expression expanding it in powers of Q. This expansion is perturbative near the critical temperature where 
Q is expected to be small. In the end one obtains the following variational expression: 

Fi2(Q,e) = i^para(/3l) + i^para(^2) + 

- nTVQ?-^TVQi-n2TVp2-^TVQ3 + 

2 2 o 

n 

- l^^^Q' + IJ:QI^QIc- ^EQt,-e^Paa + 0{Q% (8) 

abc ah a=l 

where -Fpara(/3i) and i^para(/32) are terms that do not depend on Q and are irrelevant for the present discussion. The 
above expression is 0{n) and we have neglected terms O(n^) that can be present and are relevant to evaluate the 
free energy fluctuations [24]. For the time being we assume that the only dependence on the temperature is in the 
reduced temperature ti, T2 and T12 while the other coefficients do not change with the temperature. 

The variational action in presence of a forcing term e has been computed in [25] and will be also reconsidered in 
appendix A. It turns out that the chaos effect depends crucially on the parameters w, v and C12 that is the coefficient 
of the term (n — T2)^ in the expansion of ri2: 

ri2 = \{Ti + T2) + ^(ri - T2f + O(r^). (9) 

Following [7, 25] we report the following value of the free energy shift at leading order in q and (n — T2): 

AFi2{q) = A \q\\Ti - T2f ^=^{^- C12) , (10) 

The above expression holds when q is small but it is larger than [n — T2]; in the opposite situation {q « n — T2\), 
as we will show in appendix A, we have at leading order: 

^1/2 / „ x3/2 

^FM = Bq\i-r2? B=-^{--ci2) (11) 
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A peculiar feature of the SK model is that the quantity — C12) vanishes because uj = v = C12 = 1 therefore chaos 
is not present at this order. In this case a more refined computation is necessary [5] and it shows that relationship 
(10) must be replaced with the much smaller expression 

AFM=^^\q\'AT\ (12) 

As a consequence chaos in temperature in the SK model is exceedingly weak and it was not observed in numerical 
simulations up to quite large system sizes [3]. 



IV. CHAOS IN TEMPERATURE ON BETHE LATTICE SPIN-GLASS MODELS 

A. Setting up the computation 

In [24] we have obtained the mapping of the free energy of spin-glass models defined on Bethe lattices with finite 
connectivity on the action (8). In this section we will extend those results to study chaos in temperature in these 
models. 

Extending the treatment of [26] to the case of two coupled systems we express the free energy as a variational 
functional of the order parameter p{{ai,a2}) that is a function defined on 2n Ising spins {ai} = a\,...,a^ and 
{CT2} = (72, . . . , • The variational expression of the free energy reads: 



A2(e,/3i,^2) 



M , 

— In Tr 

n 

M + 1 



2n 



exp 



In j Tr{,^,,,}Tr{,,,,,}p^{ai,a2}p^{Ti,r2} x 



/3i J ^ <Tf + /32 J ^ a2"r2" + e ^ a^a^^ + ^ ^ 



(13) 



Where M + 1 is the connectivity of the lattice and the square brackets mean average with respect to the distribution of 
J. The above expression has to be extremized with respect to p{ai,a2}- We note that it is invariant under a rescaling 
of p{a-i,a2} so that we can choose any normalization for it. If we normalize p{ai,a2} to one the corresponding 
variational equation in terms of p{a) reads: 



pWi,(^2} = jj:T^{ruT2}P {n,T2}(exp 



a=l 



(14) 



"fci 3i 



where AA is a normalization constant. 

In order to build an expansion in the order parameter we write: 

n 

^{0-1,0-2}= ^ &fei,fc2 ^ gai...afc,,ft.../3fc2^r 

ki=0,k2=0 (ai...afcj,(/3i.../3fc2) 

with 

bk„k2 = (cosh" ;3i Jcosh" ^2 Jtanh'^^ ^1 J tanh'^^ 02 J) 
The variational equation (14) can now be written as equations for qai...aki,i3i...0k2- 

Tr {.„.,}Tr ... rf^ .. r2^'^- /^{n , r2} exp e E:.i Tfrf 



9ai.. 



T: {.„r.}P^{Ti, T2} exp e ^^=1 rf 



(15) 



(16) 



(17) 



The r.h.s. of the above equations can be expanded in powers of qai...akj^,i3i...0k2- 

In general the equations for two-index qab depend on higher order objects and this fact leads to rather complex 
equations. However near the critical temperature the Q with an high number of indices terms are much smaller that 
the one with two indices and using the corresponding variational equations they can be eliminated [23]. In this way 
we can obtain an expression that depends only on a 2n x 2n matrix Q like the one appearing in the variational action 
(8). This was done in [24] for the single system and we have extended that computation to the coupled system case. 
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If we start from the equations in appendix B of [24], we divide all the equations by a factor 2 and we make a 
rescaling Qi Qi/(6ii(M - 1)), Q2 Q2/(&22(Af - 1)), P -)> P/(6i2(M - 1)), we obtain that the equations (17) 
for Q are the same that would be obtained from a variational action of the form (8) with the following coefficients: 

= 26n(M-l) (^^^ 

Mb22 - 1 

= 2622(M-1) ^''^ 
- 26i2(M-l) 

M , , 

^ = M^I (21) 
^ M{Mb, + M-2) 
(1-M64)(M- 1)2 
- _ M(M(2M-l)&4 + M-2) 

"~ (1-M64)(M-1)2 ^ "^^ 

where 64 = ((tanh/JcJ)"^) and (3c is the inverse critical temperature that obeys the equation 1 = M(tanh^ /3cJ)- The 
mixed reduced temperature fi2 is such that its expansion in terms of n and f2 is of the form (9) with the following 
expression for Ci2- 



2(M-1) (M- l)J2(l-tanh"/3eJ)2 

C12 = — — „ = (,^4j 

M (Jtanh/3cJ- Jtanh^^cJ)2M2 

The Sherrington-Kirkpatrick limit is recovered sending M to infinity and the coupling strength to zero as = 1/M. 
In this limit we have bij = l/{MTiTj), the critical temperature goes to 1 and the various coefficients read: 

n = ^r^, f2 = ^r^, fi2 = — (25) 



Cl2 = 1- 



(26) 



These are precisely the coefficients obtained for the SK model, sec e.g. [25]. The rescaling Q\ Qi/(&ii(M — 1)), 
Q2 — > Q^l^i^iM ~ 1))' P ^ Pl(b\2{.^i — 1)) was performed in order to get rid of the temperature dependence 
in all coefficients other than fi, f2 and fi2 and corresponds in the SK limit to the usual rescaling Qi — )• Qi//3i, 
Q2 — >■ Q2/ 02i P PI{PiP2)- The above definitions are such that at finite M the reduced temperature goes to 1/2 at 
zero temperature. The actual dependence of the reduced temperature with respect to the temperature is such that: 



(Jtanh/3c J - Jtanh^ 
{M - 1)T2 



(Te - T) + 0(Te - Tf 



(27) 



and the prefactor goes to 1 in the SK limit. 

The above coefficients however cannot be put simply into eqs. (10) and (11) in order to obtain the free energy 
shifts. We must bear in mind that once the variational equations (17) are expanded in powers of Q they look like as 
if they where obtained from a variational free energy of the form (8) with coefficients that in order to avoid possible 
confusion we represent as tilded. This does not means that the true free energy has an expansion with the same 
coefficients; as it was shown in appendix B of [24] this can be understood noticing that the equation for the order 
parameter corresponds to the following expression: 



= Tr 



Tr^p^(T)(exp JEe^cTc) 

Trp^(r) 



while the equation one obtains by differentiating expression (13) corresponds to: 

M-i/r l^ / / X Tr^p^(r)(expJ^^aere) 



= Tr 



({(7})aaC76 p(cr) 



Trp^(T) 



(28) 



(29) 
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Thus the two expressions arc equivalent in the sense that they have the same solution at the order at which they are 
valid. In order to obtain the expansions in powers of the order parameter matrix one could expand directly the free 
energy (13), but technically it is much simpler to expand the variational equations (17). 

This problem can be bypassed noticing that the derivatives of expressions (10) and (11) with the tilded coefficients 
allows to determine g' as a function of e. In the two different regimes we have: 

^(^-c,.j(n-r.)^ = e (30) 

and 



^1/2 /- \3/2 

In-r2|3 = e (31) 

We must take into account that the overlap q appearing in the above equation is not the true overlap. This is due to 
two reasons: 

• We must recall that we have done the rescaling Qx Qi/(6ii(M- 1)), Qi Q2/{b2i{M- 1)), P P/{h\2{M- 
1)) 



A more subtle reason is that the true overlap is given by: 

^ Tr {r,,r.}rlTlp^+'{Tu T2] exp e E^^i r^r^ 
^ {r,,r,}p''+^{ri,r2} expeE:=i rfrf 



(32) 



The difference is that there is a term while in equation (17) there is a power p^. As usual the overlap 

entering in the cavity equations is not the true overlap. 

The net effect is that in order to obtain at leading order the relationship between the true overlap and the forcing e 
we have to make the following rescaling in eqs. (30) and (31), 

6i2(M-l) 

In the SK limit the above rescaling reduce to g — > j3il32q- Near the critical temperature we have: 

fci2(M-l) M-1 



1 + 6i2 M + 1 



(34) 



The corresponding expressions yield the overlap as a function of the forcing and can be integrated back to get the 
correct free energy shifts in the two regimes (n — f2) <C g and q^ (fi — f2): 



M- 1^ ^ ~~ 



AF,2(,) = A i,|3(n -f2f ^ = ( irn J 6^ I ^ - ) ' 



and 



We can now apply the previous formulae to different distributions of the J. 



B. Diluted bimodal Distribution 



A surprising feature of the above expressions is that a direct computation shows that the quantity vjGp' — c\2 
vanishes in the case of a bimodal distribution of the coupling J = ±1 as in the SK model. Therefore for these models 
we expect chaos to be a much smaller effect possibly of the same order of the SK model. This is consistent with the 
fact that chaos is very difficult to be observed in numerical simulations of these models. 
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This can be seen considering the case of the random-bond bimodal distribution where any couphng in the lattice is 
zero with probability p or ±1 with probability 1 — p: 

P{J)=p6{J) + ^^^{S{J + l) + 5{J-l)) {0<p<l). (37) 
In this case the relevant parameters to be insert in eqs. (35) and (36) read: 

Be = arctanh — . (38) 



M _ M{1- M\l-p)-2Mp) 



M-1' (M- 1)2(1 -M(l-p)) 



(39) 



M{M{l-p) - l)arctanh' 



2 1 



In the last expression we have used eq. (27). The chaos prefactor is: 

^-^'' = M(l -p) - 1 ^^'^ 
and we see that it vanishes in the purely bimodal case corresponding to p = 0. 

C. Poissonian Distribution of the Connectivity 

The parameters computed above can be also used to study the model where the connectivity of each spin has a 
Poissonian distribution. In general we have to take the M — > oo limit in the above expressions sending p — )• 1 as 

p=I-^. (42) 

where a is the average connectivity of a site. In the case where the couplings strength is ±1 the relevant parameters 
read: 

/3c = arctanh (43) 

u = 1 (44) 
u = ^ (45) 



a- 1 
1 



^-c,2 = (46) 



|fi - fal = ^^arctanh^ -1= |Ti - T^l + 0{T^ - . (47) 



^ - ci2 = (48) 

w"' a — 1 



The free energy differences read: 



AFi2(g) = \q\' ^— arctanh^ (49) 



2 



. ^ , , 2 (2 + a)i/2(a- 1) 1 
AFi2(? = t/tt^ ^ arctanh^ ^ 



(50) 
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V. CONCLUSIONS 



We have shown that mean-field spin-glass models defined on random lattices with finite connectivity display chaos 

in temperature. Chaos is stronger than in the SK model by fom* orders of magnitude in perturbation theory. 

In general if we consider the region of small overlaps we could expect a chaos effect such as AFi2{q) oc g^|Ti — 
instead as we have seen in section (III) the effect is smaller AFi^iq) oc q'^\Ti — T^]^. In other words chaos in the 
generalized SK model, to which the random lattice models can be mapped, is larger than in SK but is nevertheless 
not as strong as one could naively expect. This is because in Bcthc lattice models chaos is truly a RSB effect [6]: a 
model with a stable RS phase has a single stable thermodynamic state that can be followed increasing or decreasing 
the temperature and therefore is not chaotic [33]. The connection between chaos and RSB is reflected by the fact 
that the coefficients of the free energy shifts cqs. (10) and (11) depend on the coefficient u of the quartic interaction 
that is responsible for RSB. 

In diluted models chaos in temperature is considerably stronger than in SK with the notable exception of Bethe 
lattices with bimodal interactions as shown by eq. (41). In this case we expect the effect to be of the same order of 
magnitude of the SK model. In the SK model one can be prove that the quadratic terms in q^ in the free energy shift 
vanishes at all orders [9, 25], we mention that the same result can be proven at all orders in the case of the Bethe lattice 
with purely bimodal interaction but this will be published elsewhere. The argument relies on the local homogeneity of 
the Bethe lattice with bimodal interactions that is also responsible of non-Gaussian free energy fluctuations [24, 32] . 

The fact that chaos in temperature on Bethe lattice with bimodal interaction is as weak as in the SK model is 
supported also by existing numerical results. The numerical data of Billoire and Marinari (BM) suggest absence of 
chaos in SK [3] , only a close look at the function Pl^^ (gr) gives a hint that the effect may be present due to a very 
slow increasing of the weight in P^'^ (0) with the system size. On the other hand the theoretical value computed in 
[5] shows that the effect is exceedingly small in the SK model and that it is practically unobservable at the system 
sizes simulated in [3] . BM considered also the Bethe lattice with connectivity c = 6 and bimodal distribution of the 
coupling finding again no strong chaos effect in agreement with the results presented here. 

Finally we would like to remark that the extension of our computations to finite dimensional models (at least for 
large dimensions) can be done using the techniques of [28]. It would be very interesting to study if we can obtain 
reliable predictions in high dimensional models, e.g. D = &. 



Appendix A: The Quadratic Action 

In this subsection we compute the free energy shift when e is much smaller than any other parameter in the theory. 
In this case the problem can be solved expanding the variational action (8) at second order around the solution 
Pab = 0. A similar treatment was recently put forward in the case of bond chaos [7, 31]. The P-dependent term in 
the action can be written as: 

Fi2(P,e) = iTrPAP-eTrP + o(p2) (Al) 

where the matrix A is given by: 

Aab = -(2ri2 -^TrQl- ^Tr Ql)6ab - w(Qi + Q2)ab -v{QI + QI + QMab (A2) 
extremizing with respect to P the above expression at given e we easily obtain: 

Pi2(e,/3i,/92) = -^TrA-i (A3) 

In order to compute the trace we diagonalize the matrix Aab- The eigenvalues of a hierarchical matrix described by 
{ad, a{x)) are given by 



J a{x 



Ao = ct(j — / a{x)dx deg : 1 

Ax = Id - [ xa{x) + / q{x)dx ] deg : - 



and the trace is given by: 



lim - y ^ = - C a{y)-^ (A4) 
n^on^Xa Al io A(y)2 ^ > 
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Let us examine the eigenvalue Aa(0): 

-A„(0) = 2ti2 - |Tr Qf - ^Tr Q\ - uj{q^ + q^) + v{ql + ql + q{q^) (A5) 

at this point we exploit the fact that when ri = T2 this eigenvalue must vanish obtaining the condition: 

Ti - ^Tr - + \vq\ = (A6) 
n 2 

summing the above equation for ri and T2 to the expression of Aa(0) we obtain 

A„(0) = -2ti2 + Ti + r2 + - q^f = -\ (ci2 - ^) (n - (A7) 

Thus we encounter the same factor (ci2 — ^) of expression (10). In order to complete the computation of the trace 
we need the expression of \a{x) around x = 0. In order to do this we use the relationship A„(x) = —xa,{:r), where 
the dot means derivative with respect to x. At leading order in ti and T2 we have a{x) = —uj{qi{x) + q2{x)). The 
solution of the free problem is such that qi{x) = q2{x) = ujx/2u [20], therefore 

2 2 

K{x) = A„(0) + {Xa{x)f = Xl{0) + A„(0)— (A8) 

2u u 

thus at leading order for small Aa(0) we have: 

lim -Tr^-^ = / 5 — ~ — — / = — — (A9) 

"^0" h u A2(0) + Aa(0)f x2 uy^\l'\<d)Jo 1 + 2/' 2u1/2a3/2(0) 

therefore 
and 

,,1/2 / w \3/2 

AFi2((Z,/3i,/32) = 9'^-^ C12-— |Ti-r2|3 (All) 
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